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Abstract 

This is a follow-up to our earlier work for the energies and the charge (vector) and matter (scalar) 
\ distributions for S-wave states in a heavy-light meson, where the heavy quark is static and the light 

quark has a mass about that of the strange quark. We study the radial distributions of higher 
angular momentum states, namely P- and D-wave states, using a "fuzzy" static quark. A new 
improvement is the use of hypercubic blocking in the time direction, which effectively constrains 
the heavy quark to move within a 2a hypercube (a is the lattice spacing). 

The calculation is carried out with dynamical fermions on a 16 3 x 32 lattice with a « 0.10 fm 
generated using the non-perturbatively improved clover action. The configurations were generated 
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by the UKQCD Collaboration using lattice action parameters /3 = 5.2, cg-^y = 2.0171 and k = 
0.1350. 

In nature the closest equivalent of this heavy-light system is the B s meson. Attempts are now 



q_) ' being made to understand these results in terms of the Dirac equation. 
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I. MOTIVATION 



There are several advantages in studying a heavy-light system on a lattice. Our meson 
is much more simple than in true QCD: one of the quarks is static with the light quark 
"orbiting" it. This makes it very beneficial for modelling. On the lattice an abundance of 
data can be produced, and we know which state we are measuring - the physical states 
can be a mixture of two or more configurations, but on the lattice this complication is 
avoided. However, our results on the heavy-light system can still be compared to the B s 
meson experimental results. 

II. MEASUREMENTS AND LATTICE PARAMETERS 

We have measured both angular and radial excitations of heavy-light mesons, and not 
just their energies but also some radial distributions. Since the heavy quark spin decouples 
from the game we may label the states as L± = L ± |, where L is the angular momentum 
and ±| is the spin of the light quark. 

The measurements were done on a 16 3 x32 lattice with dynamical clover fermions. We 
have two degenerate quark flavours with a mass that is close to the strange quark mass. The 
lattice configurations were generated by the UKQCD Collaboration. Some details about the 
different lattices used in this study can be found in Table [H Two different levels of fuzzing 
(2 and 8 iterations of conventional fuzzing) were used in the spatial directions to permit the 
extraction of the excited states. 





# of configs. m q k a [fm] m n [GeV] 


DF3 
DF4 


160 l.lm s 0.1350 0.110(6) 0.73(2) 
119 0.6m s 0.1355 0.104(5) 0.53(2) 



TABLE I: Lattice parameters. These are UKQCD Collaboration's lattices with (3 = 5.2 and 
C sw = 2.0171. 
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III. 2-POINT CORRELATION FUNCTION 



The 2-point correlation function (see Figure [T] a) is defined as 

C 2 (T) = (P t TG q (x,t + T,t)P t+T r^U Q (x,t,t + T)) , (1) 

where t/ Q (x, t,t + T) is the heavy (infinite mass)-quark propagator and G g (x, t + T,t) the 
light anti-quark propagator. P t is a linear combination of products of gauge links at time t 
along paths P and T defines the spin structure of the operator. The (...) means the average 
over the whole lattice. The energies (mj) and amplitudes (ai) are extracted by fitting the 
C*2 with a sum of exponentials, 

C 2 (T) « J2 ^~ miT , where N max = 2 - 4, T < 14. (2) 

i=i 

Fuzzing indices have been omitted for clarity. 

IV. SMEARED HEAVY QUARK 

We introduced two types of smearing in the time direction to allow the stationary quark 
to move a little, but not too far, from its fixed location. First we tried APE type smearing, 
where the original links in the time direction are replaced by a sum over the six staples that 
extend in the spatial directions (in Fig. [2] on the left). This smearing is called here "Sum6" 
for short. To smear the static quark even more we then tried hypercubic blocking, again 
only for the links in the time direction (in Fig. [2] on the right). Now the staples (the red 
ones in Fig. [2]) are not constructed of the original, single links, but from staples (the blue 




(x,t) (x, -t 2 ) 

FIG. 1: Two- and three-point correlation functions 
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FIG. 2: APE smearing in the time direction (on the left) and hypercubic blocking (on the right). 

ones in Fig. [5]). This allows the heavy quark to move within a "hypercube" (the edges of 
the "cube" are 2a in spatial directions but only one lattice spacing in the time direction). 
This smearing is called here "Hyp" for short. Smearing the heavy quark was expected to 
improve the measurements, particularly radial distributions - which it did, to some extent. 

V. ENERGY SPECTRUM AND SPIN-ORBIT SPLITTINGS 

The energy spectrum obtained is shown in Fig. [3j Using different smearing for the heavy 
quark does not seem to change the energies too much - except for the P+ state. It is not 
understood yet why this state should be more sensitive to changes in the heavy quark than 
the other states we have considered. The energy of the DH — state was expected to be near 
the spin average of the D— and D+ energies, but it turned out to be a poor estimate of this 
average. Thus it is not clear whether or not the FH — energy is near the spin average of the 
two F-wave states, as was hoped. Our earlier results can be foun in Ref. 1|. 

One interesting point to note here is that the spin-orbit splitting of the P-wave states 
is small, almost zero. We extracted the energy difference of the P+ and P— states in two 
different ways: 

1. Indirectly by simply calculating the difference using the energies given by the fits in 



4 



2- 



< DF3 

V DF3sum6 
> DF3hyp 
- fit 



1 st excited state 



i 






2S 



P- 



P+ 



D+- 



D- 



D+ 



F+- 



FIG. 3: Energy spectrum of the heavy-light meson. Here L+(— ) means that the light quark spin 
couples to angular momentum L giving the total j = L ± 1/2. 2S is the first radially excited L = 
state. The DH — is a mixture of the D— and D+ states, and likewise for the FH — . Energies are 
given with respect to the S-wave ground state (IS). Here ro = 0.525 fm was used to convert the 
energies to physical units. The error bars shown here contain statistical errors only. The solid 
line labelled as "fit" is from a model based on the one-body Dirac equation - see section IVIIII for 
details. 

Equation [2J when the P+ and P— data are fitted separately. 

2. Directly by combining the P+ and P— data (taking the ratio) and fitting everything 
in one go with 



and P+ ground states, respectively. The energy difference, rather than the energies 
themselves, were varied in the fit. The second exponential contains the corrections 
from the excited states. 

D-wave spin-orbit splitting was also extracted in a similar manner. The results are shown 





where A, B and mc are fit parameters. mp_ and mp_|_ are the energies of the P 



in Fig. HI 
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FIG. 4: The Spin-Orbit splittings of P-wave and D-wave states. The P-wave spin-orbit splitting 
seems to be small (could be almost zero), whereas the D-wave spin-orbit splitting is larger. This 
is not fully understood yet, but the numbers are still preliminary. 

VI. RADIAL DISTRIBUTIONS: 3-POINT CORRELATION FUNCTION 



For evaluating the radial distributions of the light quark a 3-point correlation function 
shown in Fig. [T] b is needed. It is denned as 



C 3 (R,T) = (^ U Q T G q i 6 G q2 ). 



(4) 



This is rather similar to the 2-point correlation function in Eq. [TJ We have now two light 
quark propagators, G q \ and G q 2, and a probe Q(R) at distance R from the static quark (74 
for the vector (charge) and 1 for the scalar (matter) distribution). 

Knowing the energies rrii and the amplitudes a; from the earlier C2 fit, the radial distri- 
butions, x y (ii!)'s, are then extracted by fitting the C 3 with 

%ax 

C 3 (R,T)^ J2 aie- mit " x ij (R)e~ m ^a j . (5) 

The results are plotted in Figures EH3 The error bars in these figures show statistical errors 
only. See [21 for earlier S-wave distribution calculations. The "Sum6" distributions have 

n 

been published in |3[, but the "Hyp" results are still preliminary. We are currently trying 
to improve the analysis of the D-wave radial distribution data. 
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VII. CHARGE SUMRULE 



When measuring a radial distribution it is easy to also measure the sumrule by simply 
summing over the whole lattice. Our results for the charge sumrule are shown in Figure [71 
We included the vertex correction Z v = 0.7731, so that with our normalisation the result 
should be one. This comes out very nicely for the S-wave, but for the other states we get 
the sumrule to be somewhat smaller. 



VIII. A MODEL BASED ON THE DIRAC EQUATION 



A simple model based on the Dirac equation is used to try to describe the lattice data. 
Since the mass of the heavy quark is infinite we have essentially a one-body problem. The 
potential in the Dirac equation has a linearly rising scalar part, bscR, as well as a vector 
part by ec R. The one gluon exchange potential, «OGE ' ^OGE> ^ s m °dified to 

_ x k 2 + Am 2 g 



V QGE (R)<x j dk Jo (kR)\n 
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FIG. 5: On the left: The S-wave ground state charge distribution. Here we compare two lattices, 
DF3 and DF4. The essential difference between the two lattices is the mass of the light quark. 
However, the effect of the light quark mass on the distribution seems to be negligible. The label 
"model" on the solid line refers to the model presented in section IVIII1 On the right: The S-wave 
ground state and 1st excited state charge distribution overlap. Note that we see one node, as 
expected from the Dirac equation. 
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FIG. 6: On the left: The P-wave ground state charge distributions. The P+ distribution has a peak 
a bit further out than the P— , which is expected. These measurements were done using a static 
heavy quark. The solid lines are predictions from the model in section IVlIII On the right: The P+ 
ground state charge distribution. Here the results with APE smeared heavy quarks are compared 
to the ones with a strictly static quark. Smearing seems to slightly improve the measurements. 
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FIG. 7: On the left: The D— ground state charge distribution. When the two different smearings 
of the heavy quark (APE type smearing, or "Sum6" for short, and hypercubic blocking, or "Hyp") 
are compared, we see that "Hyp" gives only a small improvement over the "Sum6". On the right: 
The charge sumrule. 

where Aqqj^ = 260 MeV and the dynamical gluon mass m g = 290 MeV (see [l] for details). 
The potential also has a scalar term muL(L + 1), which is needed to increase the energy of 
higher angular momentum states. However, this is only a small contribution (about 30 MeV 
for the F-wave). 

The solid lines in the radial distribution plots are predictions from the Dirac model fit 
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with m = 0.088 GeV, a QGE = 0.81, b sc = 1.14 GeV/fm, 6 V ec « 6 SC and u = 0.028. These 
are treated as free parameters with the values obtained by fitting the ground state energies 
of P-, D- and F-wave states and the energy of the first radially excited S-wave state (2S). 
Note that the excited state energies in Fig. [3] were not fitted. This fit was done using the 
energies obtained with APE smearing ("Sum6"), and the latest "Hyp" data was not used. 

IX. CONCLUSIONS 

• There is an abundance of lattice data, energies and radial distributions, available. 

• The spin-orbit splitting is small and supports the symmetry 6yec = ^sc as proposed 




• The energies and radial distributions of S-, P- and D-wave states can be qualitatively 
understood by using a one-body Dirac equation model. 

Acknowledgements 

I am grateful to my supervisor A. M. Green and to our collaborator, Professor C. Michael. 
I wish to thank the UKQCD Collaboration for providing the lattice configurations. I also 
wish to thank the Center for Scientific Computing in Espoo, Finland, for making available 
resources without which this project could not have been carried out. The EU grant HPRN- 
CT-2002-00311 Euridice and financial support from the Magnus Ehrnrooth foundation are 
also gratefully acknowledged. 



[1] UKQCD Collaboration, A. M. Green, J. Koponen, C. Michael, C. McNeile and G. Thompson, 

Phys. Rev. D 69, 094505 (2004) 
[2] UKQCD Collaboration, A. M. Green, J. Koponen, C. Michael and P. Pennanen, Phys. Rev. D 

65, 014512 (2002) and Eur. Phys. J. C 28, 79 (2003) 
[3] UKQCD Collaboration, A. M. Green, J. Ignatius, M. Jahma, J. Koponen, C. McNeile and C. 

Michael, PoS LAT2005 (2005) 205 
[4] T. A. Lahde, C. Nyfalt and D. O. Riska, Nucl. Phys. A674, 141 (2000) 



9 



[5] P. R. Page, T. Goldman and J. N. Ginocchio, Phys. Rev. Lett. 86, 204 (2001) 



10 



